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I.   INTRODUCTION 

Every  product  eventually  fails  and,  sooner  or  later,  is 
no  longer  usable  for  its  original  purpose.   The  rate  at  which 
this  deterioration  process  proceeds  is  called  the  failure  rate. 
One  is  interested  in  predicting  the  time  to  failure  for  a 
component  or  an  entire  system.   It  is  impossible  to  predict 
the  time  to  failure  precisely,  but  one  can  model  this  process, 
viewing  the  time  to  failure  as  a  random  variable  which  has 
a  probability  distribution  that  can  be  specified.   There 
exist  numerous  life  distributions  which  provide  realistic 
probabilistic  descriptions  of  the  behavior  of  components  or 
whole  systems. 

This  paper  is  concerned  only  with  the  exponential  distri- 
bution which  plays  a  fundamental  role  in  reliability.   This 
distribution  has  many  properties  which  facilitate  the  modelling 
process,  and  so  encourage  its  application  in  many  fields. 
Where  not  directly  applicable,  it  is  often  used  as  a  bound. 

According  to  the  exponential  distribution,  a  component 
has  a  constant  failure  rate.   An  equivalent  property  is  that 
the  remaining  life  of  the  component  is  independent  of  its  age: 
the  "memory less"  or  "no  aging"  or  "no  wear"  property.   This 
means  that  a  used  component  is  always  as  good  as  a  new  one. 

Convolving  life  distributions  for  components  to  get  the 
final  life  distribution  for  the  system  is  rather  difficult, 


even  for  the  exponential  distribution.   For  an  increasing 
number  of  components  it's  no  longer  possible  to  compute  a 
closed  form  life  distribution.   This  paper  will  show  how  to 
use  a  reliability  shorthand  notation  to  get  analytical  results 
for  small  basic  systems,  and  later  use  an  algorithm  to  com- 
pute numerical  results  with  an  Apple  11+  microcomputer.   The 
Pascal  programs  are  provided  in  the  appendices. 


II.   RELIABILITY  SHORTHAND 

A.   NOTATION  AND  DEFINITIONS 

The  shorthand  notation  "EXP [A,]  +  EXPtA-]"  denotes  the 
distribution  for  a  random  variable  T-.+T2.   T,  has  an  exponen- 
tial distribution  with  failure  rate  A,  and  T„  has  an  exponen- 
tial distribution  with  failure  rate  A„.   T,  and  T-  are  assumed 
to  be  independent. 

In  the  following  T  always  denotes  a  random  life  length 
and 

F(t)   =   Probability  (T  >  t)   =   1  -  F(t) 

is  the  survival  function  as  a  function  of  time. 

The  convolution  of  life  distributions  corresponds  to 
the  summation  of  components  with  independent  random  lives 
to  get  the  life  distribution  of  a  system. 
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III.   RELIABILITY  FOR  A  ONE  COMPONENT  SYSTEM  WITH  SPARES 

A.   GENERAL  DESCRIPTION 

There  are  two  ways  for  a  simple  system,  consisting  of 
one  active  component  and  one  spare,  to  survive  a  specified 
time  t.   Either  the  active  component  survives  time  t,  in 
which  case  the  spare  need  not  be  used,  or  the  active  com- 
ponent fails  at  some  time  s  (0  <_  s  <_  t)  and  the  spare  takes 
over  and  doesn't  fail  in  the  interval  from  s  to  t  and  thus 
completes  the  mission. 

Let  F,  (t)  and  F_(t)  be  the  exponential  survival  functions 
for  component  one  and  two  respectively,  and  f,(t)  and  f?^) 
their  respective  densities.   Then  the  survival  function  F(t) 
for  the  system  is  denoted  by: 

t 
F(t)   =   F  (t)  +   /   F  (t-s)f, (s)ds 
L  0  Z 

where  the  integrand  shows  the  probability  that  component  one 
fails  at  time  s  (f, (s))  and  component  two  takes  over  for 
at  least  the  length  of  time  t-s.   The  following  shows  some 
examples  how  to  get  the  survival  function  in  closed  form  for 
different  basic  convolutions. 
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B.   BASIC  SURVIVAL  FUNCTIONS 

1 .   System  with  N  Identical  Components 

The  shorthand  notation  for  this  system  is: 

EXP  [A]  +  EXP  [A]  +  ...  +  EXP  [A  J 

Whenever  an  active  component  fails  a  spare  takes  over,  as 
long  as  there  is  a  spare.   When  the  n-th  component  fails,  the 
whole  system  stops  functioning. 

For  n  =  2 : 


t 
F(t)   =   F,  (t)  +   /   F  (t-s)f,  (s)d: 

■  0 


e"Xt  4-   /*  e-A(t-s)Ae-Asds 
0 


(l  +  At)e"At/  t    0 


The  density  function  for  this  survival  function  is 
f(t)   .   .  dFl(t)   _     d[(l+At)e-At] 


dt  dt 


-[(1+At) (-Ae~At)  +  e"AtA] 


,  ,  -At   ,2.  -At  ,  ,  -At, 
-  [-Ae    -  A  te    +  Ae 


X2te~Xt 
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For  n  =  3 : 

The  survival  function  for  the  two  component  system 
will  be  convolved  with  the  survival  function  of  the  next 
component  to  obtain  the  survival  function  for  the  three 
component  system. 


F(t)   =   (l+At)e'At  +   /   e"A(t's)A2se'Xsds 

0 


which  leads  to: 


2 

F(t)   =   (1  +  At  +  (A2t,)  )e~At,        t  >  0. 


For  all  n: 

Continuing  in  this  manner  one  gets  the  general  solution 
for  the  survival  function  of  this  kind  of  system: 


n      .  , 
I     (At)1"1 

F(t)   =   — ITE —  ,         t  >  0. 

(i-l)e 


This  class  of  convolutions  always  gives  the  well-known 
Erlang[n,A]  survival  function. 

2 .   System  with  N  Dissimilar  Components 

The  shorthand  notation  for  this  system  is  denoted 
by: 


EXP [A,]  +  EXP[A2]  +  ...  +  EXP[AnJ 
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For  n  =  2 : 

The  survival  function  for  the  active  component  is 


-A,t 
F,(t)   =   e      ,  t  >  0 


The  survival  function  for  the  spare  is 


-At 
F2(t)  =      e      Z       ,  t    >_   0 


The   survival   function   for   this   system   is 


t 
F(t)       =      Fx(t)    +      /      F2(t-s)f1(s)ds 


-A,t  t   -A^(t-s)         -A.s 

1< 


e1      +      /e2  A„e     1  ds 


0 


-A,t  A,  -At         t  -(A    -A    )s 

e      x      +  t — ^— e      z        j    (A,-A_)e        x      z      ds 
Al  "  A2  0  l 

-At              X  -At  -<VA2)t 

e  +  t r—  e  (1   -   e  ) 

Al"  A2 

—  A  .  t  — At  —  A  .  t 

(A1-A2)e  +    X±e      2      -    X^  ) 

uq  -  a2) 

-At  -At 

(A,e      Z      -    A0e  ) 

1 t t    >    0 

(A,    -    X)  ■  t-°- 
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For  n  =  3 : 

The  resulting  survival  function  for  the  two  component 
system  will  be  convolved  with  the  survival  function  of  the 
next  component  to  obtain  the  survival  function  for  the  three 
component  system. 

The  result  is 


F(t)  A2A3        "Alt       A1A3         "A2t 

(A3-A1)  (X2-X1)  e      +  (A3-A2)  (Al"A2)  S 


A1A2        "A3t 
(X3-X2)(X  -A  )  e 


For  all  n: 

Continuing  in  this  manner  one  gets  the  general  solution 
for  the  survival  function  of  this  kind  of  system: 

n  ...    X .       -A  .  t 
F(t)   =    7   J2fi 3 e   1 

*  {   '  .s  n  (a.-a.)  e 

1=1  ■_,■      1      1 

The  two  general  survival  functions  for  systems  built  from 
components  which  all  have  dissimilar  failure  rates  or  all 
have  the  same  failure  rate  are  very  useful  for  computations 
of  survival  functions  of  systems  which  are  formed  from  sub- 
systems of  the  above  two  kinds. 

It  will  be  shown  in  the  following  section  how  to  use 
the  above  results  to  get  a  closed  form  for  the  survival 
function  of  a  system  which  can  not  be  computed  with  only  one 
of  those  general  formulas  . 
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3 .   Systems  With  Identical  and  Dissimilar  Components 
Shorthand  notation: 

EXP[A1]  +  EXP[X1]  +  EXP[A2]  +  EXP[A2J 

Let  the  random  variables  T  be  distributed  in  the  following 
way: 


T±      ~   EXP[A1],  T2   ~   EXP [A,], 


T3   -   EXP[A2],  T4   -   EXP[A2J. 


There  are  different  approaches  available  for  computing  the 
survival  function  for  this  system. 

One  can  compute  the  survival  function  for  the  first 
two  and  then  the  second  two  components,  using  the  general 
formula  of  Section  III.B.l  and  then  convolve  those  two  sur- 
vival functions;  or  one  can  compute  the  survival  function  of 
the  first  two  components  and  then  add  the  other  components 
one  at  a  time.   The  latter  approach  will  be  used  here  to 
give  a  general  idea  of  the  procedure,  and  to  lead  to  an 
algorithm  for  the  numerical  evaluation  of  the  reliability  of 
systems  consisting  of  so  many  components  that  a  closed  form 
is  no  longer  attainable  or  not  attainable  in  an  efficient 
manner. 

The  first  two  components  (T'  =  T,,  T*  =  T,  +TO 
have  as  shown  in  III.B.l  the  survival  function 
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_        _        t 

FT,(t)   =   F  ,(t)  +   /   F   (t-sjf.,,  (s)ds 
2  Ll  0     T2       Tl 

~A  ,  t 

=   (ai;L  +  a12t)e  t  >_  0  (1) 

with: 


an  =  1  and     ai2  =  V 

Adding  the  next  component  leads  to  T'  =  T,  +  T      +   TV 
and  the  survival  function 


FT'(t)   =   F   (t)  +   /   F  ,(t-s)fT,(s)ds 
2       0     l3  L2 


-A,t      t  -A  (t-s)     -A-S 

2, 

1' 


-A,  t         U 

=   (1  +  Axt)e     +   /  e   2      A?se 


0 


which  leads  after  the  integration  to 

{X22  -  2A1A2)  +  (A1A2  -  A^A0)t      -A-t  A^  -At 

F  ,  (t)   =   » e      +  ~  e 

3  (Ax-A2)z  (A1-A2)z 

or  in  the  form  of  equation  (1)  to 


-A,t       -A2t 
FT,(t)   =   (a11+a12t)e      +  a^e 


with 


2  2   2  2 

A^~4A-]  A^  A  -.  A  ~  —  A-i  A^  A  -i 

al  1   =           2'  al?   =            2~'  a?l              ~~2 

11        u1-\2)z  lz         (ArA2}  (a1-a2)z 
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Adding  the  last  component  leads  to  T '  =  T,  +  T  +  T  +  T   and 
the  survival  function: 


t 


F  , (t)   =   F   (t)  +   /   F_,(t-s)f_,(s)ds 
4  x3       0     X4       L3 


-A1t        -At 
(a11+a12t)e      +  a21e 


t  -A2(t-s)  -As 

+  J    e         [<al2-allXl-a12Xls)e 

-A„s 
-  A2a2]_e     ]ds 


Integrating  and  collecting  terms  gives  the  final  survival 
function  for  the  whole  system: 


-A,t  -At 

FT,(t)  =  (a11+a12t)e   x   +  (a21+a22t)e     ,    t  >  0 

with 


3     2  2 

A2  -  3A1A2  \1k2 


all  3  '  a12  ' 

(A2-A1)J  1Z      (A2-A1)Z 


3     2  2 

Xl  "  3X1A2  A1A2 

a01   -  =■  ,  a. 


21  3  '  22  2  * 

(A2-A1)J  AA  (AX-A2)Z 

This  section  demonstrates  that  one  can  get  the  survival 
function  by  adding  one  random  variable  at  a  time  and  thus 
computing  the  new  survival  function.   Each  summation  then 
gives  a  new  or  updated  a. ..   Even  though  the  integration  is 
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not  shown  explicitly,  one  can  imagine  that  this  method  of 
convolving  survival  functions  is  very  tedious  and  it  is  very 
easy  to  make  mistakes.   So  the  next  section  introduces  an 
algorithm  to  compute  all  the  a. .,  which  then  makes  it  easy 
to  evaluate  the  survival  function  at  any  time. 

C.   ALGORITHM  TO  COMPUTE  THE  CONVOLUTION  OF  EXPONENTIAL  LIVES 
The  shorthand  notation  to  be  solved  is 


EXP[A,]  +  EXP[X11  +  ...  +  EXP[A,]  (nx  terms) 

+  EXP[A2]  +  EXP[A2]  +  ...  +  EXP[A2]  (n2  terms) 


+  EXP[AR]  +  EXP[Ak]  +  ...  +  EXP[Ak]  (nfc  terms) 

The  general  solution  for  the  survival  function  is  of  the  form 


-A,t         -At  -At 

F(t)   =   A1(t)e  x      +   A2(t)e   z   +  ...  +  Ak(t)e     ,   t  >  0 


with 


2       .3  Jnl-1] 


Ax(t)   =  au  +  a12t  +  a13t  +  a14t  +  . . .  +  a^t 

2       3  (n2_1) 

A2(t)   =   a21  +  a22t  +  a23t   +  a24t   +  ...  +  a^t 


2         3  (nk_1) 

Ak(t)   -   akl  +  ak2fc  +  ak3t   +  ak4fc   +  '••  +  aknkfc 

An  example  shows  how  to  get  the  survival  function  for  the 
following  shorthand  notation: 
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EXP[A1]    +    EXP  [A.,]    +    EXP [A,] 
+    EXP[A    ]     +    EXP[A2] 
+    EXP[A3] 


Survival   Function: 


o       "*A.t  -  A-t  —  A_t 

F(t)       =       ^an+a12t+ai3t    )«  +    (a2i+a22t)e  +    a31e 

t    >    0. 


An  algorithm  to  compute  all  the  factors  a,,,  a,-,  .  ..,  a->, 
and  thus  make  it  possible  to  convolve  all  the  different 
exponential  functions  was  earlier  introduced  by  Gursel  [Ref. 
1] ,  and  is  given  here  in  a  slightly  modified  form.   This 
algorithm  adds  one  exponential  random  variable  in  each  run. 
The  above  example  can  be  solved  in  this  way  by  looping  six 
times  through  this  algorithm. 
Notation  for  the  algorithm: 


K   =   number  of  dissimilar  failure  rates  (3  in  above 
example) ; 

1.   =   i-th  type  failure  rate,  i  =  l,...,k; 

nn .   =   number  of  random  variables  having  the  9-th 
type  failure  rate; 

1.    =   failure  rate  of  the  ie-th  lambda,  ie  =  1,...,K 
le 

n.   =   current  number  of  random  variables  having  the 
1     i-th  failure  rate,  n.  =  l,...,nn.; 

a.  .   =   j-th  coefficient  of  the  i-th  polynomial 

ie  =   index  for  current  polynomial  and  currently 
entering  life. 
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Algorithm  to  compute  the  coefficients: 


Initialization:   a.   =  0,  for  all  j ,k  where  j  =  1,...,K, 

-*       k  =  1 ,  .  .  .  ,  nn  . 

1 

n.  =  0,  for  all  i  where  i  =  1,...,K 

Input:  A.  for  all  i  where  i  =  1,...,K 

nn .  for  all  i  where  i  =  1,...,K. 
l 

After  the  first  run:   a,,  =  1,  n,  =  1 

n.   =  n.   +1,  until  n.   =  nn . 
le    le  le     le 


1.   Update  the  coefficients  a.   ,  for  k  =  2,...,n. 

16, k  ie 


a.   ,   =   A.  *a.     ,/(m-l),   where  m  =  n.  -j  for 
ie,k      ie   ie,m-l'  ie  J 

j  =  0 , . . . ,n.  -1 

J  ie 


2.   Update  the  coefficients  a.   , 
r  ie,l 


ni   a., A .  ni  (j-1) IX.  a. . 

V    r  ll  1  V     J        ie  IT-, 

ie,l      ie,l    .n.   (A. -A.  -L0   ,,   ,   »j 

l^ie    i   ie'  D=2   (A.-A.  )J 

i=l  i^ie    x   ie 


3.   Update  the  other  coefficients  a. .  for  all  i,k  where 

13    i  <>  ie,   n.  >  0 

for  i  =  1, . . . ,K, 

k  =  1,  . . . ,n . 


a.      =   a.    A.  /(A.  -A.)     for  all  i  ^  ie, 
i,n .       i/n.  ie'   ie   i 
i         l 


a.      =   (a.    A.  -ma.   )/(A.  -A.)    for  all  i  where 
"i .        i  ,n  . 
i  i 


i,n4       i,n,  ie   i,m    ie   i     -  f   ie  and  n   >  x 


.  i 
for  m  =  n.-j , 


j  =  l,...^(ni-l) 
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D.   PROGRAM  FOR  NUMERICAL  EVALUATION  OF  THE  RELIABILITY 
FOR  A  ONE  COMPONENT  SYSTEM  WITH  SPARES 

The  above  algorithm  is  used  to  compute  the  reliability  of 

a  one  component  system  with  similar  and/or  dissimilar  spares 

for  any  given  time.   The  input  to  the  program,  which  is 

written  in  PASCAL  and  listed  in  Appendix  A,  will  be  explained 

here  with  the  same  example  as  above: 

EXP[AXJ  +  EXPUJ  +  EXP [A,] 

+  EXP[A2]  +  EXP[X2] 

+  EXP[A3] 

Let's  assume  that  A,  =  0.1,  A  =   0.2,  X.    =  0.3  and  we  want 

to  know  the  reliability  for  t  =  0,  10,  20,  ...,90  hours.   The 

prompts  of  the  program  are  typed  in  large  letters. 

"HOW  MANY  DISSIMILAR  LAMBDAS  DO  YOU  WANT  TO  USE?" 

3 

"WRITE  DOWN  ALL  THE  DISSIMILAR  LAMBDAS  WITH  AT  LEAST  ONE 

BLANK  IN  BETWEEN  (INPUT  NUMBERS  <  1  WITH  A  LEADING  ZERO)   :" 

0.1  0.2  0.3 

"WRITE  DOWN  THE  NUMBER  OF  EACH  OF  THE  LAMBDAS  WITH  AT  LEAST 

ONE  BLANK  IN  BETWEEN:" 

3  2  1 

"HOW  MANY  TIME  POINTS  DO  YOU  WANT  TO  BE  COMPUTED?  (<=20)  :" 

10 

"AT  WHAT  TIMES  DO  YOU  WANT  THE  SURVIVAL  FUNCTION  TO  BE 

EVALUATED?   (TIME  SHOULD  HAVE  THE  SAME  DIMENSION  AS  THE 

LAMBDAS  AND  SHOULD  BE  AT  MOST  20  INTEGERS  BETWEEN  0  AND  70. 
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THE  UPPER  LIMIT  70  IS  REQUIRED  ONLY,  IF  A  GRAPHICAL 
PRESENTATION  IS  DESIRED  LATER)   :" 

0  10  20  30  40  50  60  70  80  90 

"DO  YOU  WANT  A  SIMULATION  TO  BE  DONE  PARALLEL  TO  THE 
COMPUTATIONS?   (SERIOUS  SIMULATION  NEEDS  A  LOT  OF  TIME!)   HIT 

1  FOR  "YES"  AND  "0"  FOR  NO:" 
1 

"HOW  MANY  NUMBERS  DO  YOU  WANT  TO  USE  AS  BASIS  FOR  THE 

SIMULATION?  :" 

1000 

After  the  computations  the  program  will  produce  the  following 
output : 

FOLLOWING  LAMBDAS  WERE  ENTERED:   0.1  0.2  0.3 
#  OF  EACH  LAMBDA:     3    2    1 


Results: 

RELIABILITY 

!  T 

> 

0  ) 

= 

1 

.000000 

(SIM 

= 

1.000) 

RELIABILITY 

:  t 

> 

10  ) 

= 

0 

.995799 

(SIM 

= 

0.996) 

RELIABILITY 

[  T 

> 

20  ) 

= 

0 

.925259 

(SIM 

= 

0.925) 

RELIABILITY 

!  T 

> 

30  ) 

= 

0 

.739626 

(SIM 

= 

0.728) 

RELIABILITY  1 

;  t 

> 

40  ) 

= 

0 

.508910 

(SIM 

= 

0.489) 

RELIABILITY  1 

:  t 

> 

50  ) 

= 

0 

.311271 

(SIM 

= 

0.334) 

RELIABILITY  1 

.  T 

> 

60  ) 

= 

0 

.174439 

(SIM 

= 

0.173) 

RELIABILITY 

!  T 

> 

70  ) 

= 

0 

.091597 

(SIM 

= 

0.068) 

RELIABILITY 

:  t 

> 

80  ) 

= 

0 

.045784 

(SIM 

= 

0.036) 

RELIABILITY 

[  T 

> 

90  ) 

= 

0 

.022028 

(SIM 

= 

0.020) 

"DO  YOU  ALSO  WANT  A  GRPAHICAL  PRESENTATION  OF  THE  DATA?   (HIT 
"1"  FOR  "YES"  AND  "0"  FOR  "NO"):" 
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The  graphical  output  will  not  be  shown  here,  but  will  be 
generated  from  the  program  if  desired. 
Some  notes  on  the  program: 

1.  The  program  is  set  up  to  allow  20  dissimilar  lambdas. 

2.  It  is  necessary  to  type  leading  zeros  in  PASCAL  if  one 
inputs  decimal  numbers. 

3.  The  simulation  generates  n  numbers  (in  this  example 
n  =  1000)  for  every  random  variable,  to  produce  n 
observations  on  the  whole  system. 

4.  Because  of 

the  behavior  of  the  exponential  function, 

the  sequence  of  the  lambdas  in  the  input, 

the  method  of  adding  always  only  one  random  variable 

at  a  time, 

the  possible  closeness  of  the  values  of  the  lambdas 

and 

the  possibly  large  number  of  the  dissimilar  lambdas, 
the  program  might  have  to  compute  with  very  large  and 
small  numbers  and  could  produce  erroneous  results. 
Therefore  a  simulation  is  provided  within  the  program 
to  see  if  the  computations  are  close  to  the  simulation. 
Fixed  limits  which  avoid  errors  can't  be  given,  because 
all  the  possible  errors  are  closely  related.   Fewer 
dissimilar  lambdas  for  example  would  allow  that  their 
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values  could  be  closer  together,  or  that  less  care  is 
needed  with  the  sequence  in  the  input,  etc. 
But  after  a  long  test  phase  and  a  comparison  with  results, 
produced  with  the  IBM  3033,  some  hints  can  be  given: 

1.  The  input  sequence  of  the  lambdas  should  go  from  the 
smallest  to  the  largest, 

2.  There  should  be  no  more  than  20  lambdas  altogether, 

3.  The  difference  in  value  of  the  dissimilar  lambdas  should 
be  larger  than  0.05. 
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IV.   RELIABILITY  FOR  A  SYSTEM  OF  ACTIVE  COMPONENTS  AND  SPARES 

A.   MIXTURE  OF  CONVOLUTIONS 

Before  solving  systems  with  active  components  and  spares, 
the  mixture  of  convolutions  will  be  explained  briefly.   Let's 
assume  that  we  want  to  compute  the  reliability  of  a  parallel 
system  with  two  different  components,  like  that  in  Figure  1. 


Figure  1.   Simple  Parallel  System 

The  system  starts  with  a  life  distributed  as  EXPU-^A  ].   Then 
there  are  two  different  possibilities:   either  component  #1 
fails  with  probability  P,  =  A,/(A,+A2)  and  the  residual  life 
is  distributed  as  EXP [A-].   This  is  the  upper  branch  of 
Figure  2.   Or  component  #2  fails  with  probability 
P   =  A2/(A1+A2)  and  the  residual  life  is  distributed  as 
EXP[AX]  . 
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EXP[A1+A2]  + 


EXP  [  A 


A1+A2 


EXPtA-J 


Figure  2.   Shorthand  Notation  for  Parallel  System 

Using  the  "MIX"-notation,  the  life  can  be  written  in  the 
shorthand  notation  as: 


EXP[AX+A2]  +  MIX{P1*EXP[A2] ,P2*EXP[A1] }. 


Because  of  the  distributive  law,  which  is  valid  for  these 
mixtures,  one  can  transform  this  equation  to 


MIX{P1*EXP[A2]  +  EXP[A1+A2] ) ,P2*(EXP[A1]  +  EXPU-j+Aj])} 


The  first  and  the  second  expressions  in  parentheses  multiplied 
by  their  respective  probabilities,  is  a  mixture  of  convolu- 
tions.  The  following  section  demonstrates  how  those  mixtures 
can  be  solved  numerically  even  if  an  analytical  solution  is 
not  feasible. 
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B.   DESCRIPTION  OF  THE  SYSTEM 

This  section  is  concerned  with  computing  the  reliability 
of  a  system  of  components.   The  different  components  might 
be  supported  by  spares,  which  have  the  same  or  a  different 
failure  rate  as  the  original  component.   The  following  example 
will  be  used  to  illustrate  the  idea  of  describing  the  system 
in  the  reliability  shorthand  notation.   The  listing  of  a 
PASCAL  program  for  these  types  of  systems  is  provided  in 
Appendix  B.   The  program  incorporates  the  method  of  program 
#1  to  add  the  components,  one  at  a  time,  to  get  first  the 
convolution  of  the  exponential  distributions  and  then  via 
the  mixture  the  reliability  of  the  whole  system. 
Example: 

The  system  consists  of  two  parallel  components  (#1  and 
#2)  in  series  with  a  third  component  (#3) .   A  spare  is  pro- 
vided for  component  #1.   The  components  have  the  failure  rate 
A,,  \~    and  X^   respectively.   The  spare  has  the  same  failure 
rate  as  the  component  it  can  replace. 


/• 


Figure  3.   Combined  Serial/Parallel  System  with  Spare 
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C.   SHORTHAND  BRANCHING  DIAGRAM 


x         /Exp[x1+x2+x3]- 


A1+A2+A3 


h+X2+h]+tx-^xf^[h+X3]  + 


A1+A2+X3 


\ 


X+X+A  /^W 
/ 

^^-^EXPtX^j 
ZEBO 


A1+A2+A3 


j-^-^EXPtX  +i  ] 


ZERO 


(1) 


(2) 


(3) 


(4) 


(5) 


ZERO   (6) 


Figure  4 .   Shorthand  Branching  Diagram  for  Serial/ 
Parallel  System  with  Spare 


The  numbers  in  the  numbered  parentheses  in  the  figure 
above  describe  the  different  endpoints  (leaves)  in  this  nota- 
tion.  The  connections  between  the  original  system  (root) 
and  the  leaves  will  be  called  branches.   In  the  following: 
the  connection  from  the  root  to  leaf  #1  is  called  branch  #1, 
from  the  root  to  leaf  #2  branch  #2,  etc. 

The  original  system  has  a  life  distributed  as  EXP [A,+A2+A-J 
Three  different  things  could  happen  now.   Either  component 
#1  fails  with  probability  A,/(A,+A  +A-J  or  component  #2  fails 


29 


with  probability  A2/U,+A2+A3)  or  component  #3  fails  with 
probability  A  _/  (A-,+A2+A3)  .   Let's  continue  in  branch  #1  and 
assume  component  #1  failed  with  the  above  probability.   Now 
the  spare  can  replace  the  failed  component  and  the  system's 
life  is  again  distributed  as  EXP [A , +A2+A -] .   The  system  is  now 
in  the  same  state  as  at  the  beginning  and  the  same  arguments 
hold.   Let's  assume  new  component  #1  fails  with  probability 
A ,/  (A,+A2+A  -.)  ,  then  the  remaining  life  is  distributed  as 
EXP[A2+A.J.   If  component  #2  fails  (the  probability  of  this 
failure  is  A~ (A, +A~+A 3) )  ,  the  remaining  life  is  distributed 
as  EXP[A,+A^].   The  third  possibility,  component  #3  fails  with 
probability  A  ./  ( A-,  +A2+A  _)  ,  brings  the  failure  of  the  whole 
system,  denoted  by  the  "ZERO"-notation.   This  says  that  there 
is  a  remaining  life  of  zero.   The  reason  why  this  third  branch 
isn't  simply  omitted  is  that  even  if  there  is  no  remaining 
life,  there  is  a  specified  probability  that  this  event  occurs. 
And  this  has  to  be  taken  into  account  to  compute  the  total 
life  of  the  system. 

The  other  three  branches  of  the  above  figure  are  developed 
analogously. 

D.   THE  SURVIVAL  FUNCTION 

Let's  denote  the  survival  function  of  branch  #1  by  F-^t), 
and  the  probability  that  the  events  in  branch  #1  occur  in 
this  special  sequence  as  P, .   Then  the  life  of  branch  #1  is 
distributed  as  P*  (EXP  [Aj+Aj+A-]  +  EXP  [A1+A2+A3J  +  EXPIA2+A3J  ) 
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and  the  survival  function  of  this  system  is  the  following 
mixture: 

F(t)   =   MIX{p1*F1(t) rP2*F2(t) ,p3*F3(t) ,p4*F"4(t) ,p5*F"5(t)  fp6*Fg(t) } 

For  the  computation  of  the  coefficients  for  every  single 
survival  function  the  algorithm  of  Section  III.C  can  be  used 
again. 

E.   PROGRAM  FOR  THE  NUMERICAL  EVALUATION  OF  THE  RELIABILITY 
FOR  SYSTEMS  WITH  SEVERAL  ACTIVE  COMPONENTS  AND  SPARES 

This  section  describes  the  program  for  computing  the 

reliability  for  systems  as  described  above.   While  reading 

this  part,  one  should  look  at  Figure  4.   The  program  is 

interactive.   Only  the  input  of  Figure  4  will  be  explained 

here . 


Prompt 


Input 


"How  many  branches?" 

1. 

Branch  : 

1. 

Element  : 

2. 

Element  : 

3. 

Element  : 

4. 

Element  : 

5. 

Element  : 

6. 

Element  : 

2. 

Branch  : 

1. 

Element  : 

2. 

Element  : 

3. 

Element  : 

6  CR 


(Carriage  return) 


2  CR  3  CR  0  CR 

1  CR  1  CR  2  CR  3  CR  0  CR 

1  CR  2  CR  3  CR  0  CR 

1  CR  1  CR  2  CR  3  CR  0  CR 

1  CR  2  CR  3  CR  0  CR 

0  CR 


1  CR  3  CR  0  CR 

2  CR  1  CR  2  CR  3  CR  0  CR 
1  CR  2  CR  3  CR  0  CR 
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/ 


4. 

Element 

5. 

Element 

6. 

Element 

3. 

Branch 

1. 

Element 

2. 

Element 

3. 

Element 

4. 

Element 

5. 

Element 

6. 

Element 

4. 

Branch 

1. 

Element 

2. 

Element 

3. 

Element 

4. 

Element 

5. 

Element 

6. 

Element 

5. 

Branch 

1. 

Element 

2. 

Element 

3. 

Element 

4. 

Element 

5. 

Element 

6. 

Element 

6. 

Branch 

1. 

Element 

2. 

Element 

3. 

Element 

4. 

Element 

1  CR  1  CR  2  CR  3  CR  0  CR 
1  CR  2  CR  3  CR  0  CR 
0  CR 


0  CR 

3  CR  1  CR  2  CR  3  CR  0  CR 

1  CR  2  CR  3  CR  0  CR 

1  CR  1  CR  2  CR  3  CR  0  CR 

1  CR  2  CR  3  CR  0  CR 

0  CR 


1  CR  3  CR  0  CR 

1  CR  1  CR  3  CR  0  CR 

1  CR  3  CR  0  CR 

2  CR  1  CR  2  CR  3  CR  0  CR 
1  CR  2  CR  3  CR  0  CR 

0  CR 


0  CR 

3  CR  1  CR  3  CR  0  CR 

1  CR  3  CR  0  CR 

2  CR  1  CR  2  CR  3  CR  0  CR 
1  CR  2  CR  3  CR  0  CR 

0  CR 


0  CR 

3  CR  1  CR  2  CR  3  CR  0  CR 

1  CR  2  CR  3  CR  0  CR 
0  CR 
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The  branches  in  the  figure  count  from  top  to  bottom  and  the 
elements  in  each  branch  from  the  leaf  to  the  root.   The  first 
element  in  the  first  branch  is  EXP[A_  +  A_J  and  one  inputs  only 
the  two  indices  of  the  lambdas.   The  second  element  in  the 
same  branch  is  A,/ ( A, +  A2+A-J  and  one  inputs  first  the  index 
of  the  lambda  in  the  numerator  and  then  the  indices  of  the 
denominator.   After  each  element  input  a  "0".   If  there  is  no 
sixth  element  as  in  branch  #1  to  branch  #5,  no  fourth  element 
as  in  branch  #6,  or  no  element  as  the  ZERO  elements  in 
branches  #3,  #5,  and  #6,  simply  input  a  "0".   Because  of  the 
data  structure  used  in  the  PASCAL  program,  the  system  always 
has  to  be  split  up  into  pieces,  so  that  there  is  only  a  single 
failure  in  each  step.   Because  of  the  algorithm  used,  there 
shouldn't  be  more  than  2  0  elements  in  each  branch  to  avoid 
errors. 
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V.   SUMMARY 

The  reliability  shorthand  is  a  convenient  way  to  describe 
the  life  of  different  systems.   This  shorthand  then  is  the 
underlying  data  structure  to  compute  numerically  the  relia- 
bility for  systems  for  which  it  is  difficult  or  impossible 
to  get  the  survival  function  in  closed  form. 
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APPENDIX  A 


PROGRAM 

CONVOLUTION; 

USES  TRANSCEND. 

.APPLESTUFF; 

CONST 

ASTERISKS  *'  . 

BLANK 

_  >      > 

DASH 

as*  — ■  *  > 

LINE 

—  f    i  »  . 

TYPE 

ARR  = 

PACKED 

ARRAY  CI. 

.203 

OF 

REAL; 

ARI  = 

PACKED 

ARRAYC1. 

.203 

OF 

INTEGER; 

VAR 

A: PACKED  ARRAY C 1 . . 20, 1 . . 20  3  OF  REAL; 

L,PRO,REL,RELSIM  :ARR; 

NINIT,NI,T:ARI; 

KK, JJ, IE, K, TMAX, SIMMAX, Y, X, COUNT, VAL: INTEGER; 

I, II, J,FACT,NN,NKK,NNI,CHGR,CHSIM: INTEGER; 

IT,D  .-INTEGER;    (*  INDEX  FOR  T-ARRAY  *) 

SUM, RELI: REAL; 

PR I NT VAR: CHAR; 

<*  L  =  ARRAY  OF  DISSIMILAR  FAILURE  RATES  t) 

(*  NIC  1 3  =  NUMBER  OF  LIVES  HAVING  THE  ITH  FAILURE  RATE  *) 

<*  ACI,J3  =  COEFFICIENT  FOR  THE  POLYNOMIAL  *> 

it    T  =  ARRAY  OF  TIMES  AT  WHICH  THE  RELIABILITY  WILL  BE 

COMPUTED  *> 

PROCEDURE  ATTENTION; 
VAR 

DURATION, Is  INTEGER; 


PITCHrARRAYCl. 

.  103 

OF 

INTEGER; 

BEGIN 

PITCHC13: 

=  15; 

PITCHC23: 

=  17; 

PITCHC33: 

=  19; 

PITCHC43: 

=  20; 

PITCHE53: 

DURATION: 

=  30; 

FOR  I:=l 

TO  5 

DO 

NOTE(PITCHC: 

[3, DURATION) ; 

END; 

PROCEDURE  DECIDE2; 
BEGIN 

WRITELN;WRITELN; 

WRITELNCDO  YOU  WANT  TO  CONTINUE  WITH  THE  SAME  PROGRAM? 

TYPE  "U"  : ' ) ; 

WRITELNCDO  YOU  WANT  TO  COMPUTE  MIXED  CONVOLUTIONS? 

TYPE  "XP2"  : *   ) ; 

EXIT(Pl) ; 
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END; 

PROCEDURE  INPUT; 
BEGIN 

WRITELN;WRITELN; 

GOTOXY(0, 12)  ; 

WRITELN  CHOW  MANY  DISSIMILAR  LAMEDAS  DO  YOU  WANT  TO  USE?   '); 

READLN(K) ; 

WRITELN; 

WRITELN ('WRITE  DOWN  ALL  THE  DISSIMILAR  LAMBDAS  WITH  AT 

LEAST  ONE  BLANK  IN  BETWEEN:'); 

WRITELN (' (NUMBERS  <  1  NEED  A  LEADING  ZERO!)'); 

FOR  I:=l  TO  K  DO 

READ (LCI  3) ; 
WRITELN; 

WRITE ('WRITE  DOWN  THE  NUMBER  OF  EACH  OF  THE  LAMBDAS  WITH  AT 
LEAST  ONE  BLANK  IN'); 
WRITELN ('BETWEEN: ' ) ; 
FOR  I:=l  TO  K  DO 

READ(NINITCI3) ; 
WRITELN; 

WRITE  CHOW  MANY  TIME  POINTS  DO  YOU  WANT  TO  BE  COMPUTED"1 
(  <=  20  )  :  '  )  ; 
READLN(TMAX) ; 

WRITELN CAT  WHAT  TIME  DO  YOU  WANT  THE  SURVIVAL  FUNCTION  TO 
BE  EVALUATED?' ) ; 

WRITELN C  (INTEGER  TIMES  BETWEEN  0-60  IN  THE  SAME  UNITS  AS 
THE  LAMBDAS) ' ) ; 
FOR  I s-1  TO  TMAX  DO 

READ(TCI3) ; 
WRITELN; 

WRITELN ('DO  YOU  WANT  A  SIMULATION  TO  BE  DONE  PARALLEL  TO 
THE  COMPUTATIONS?'); 

WR I TELNC  (SERIOUS  SIMULATION  NEEDS  A  LOT  OF  TIME!)'); 
WRITE  CHIT  "l"  FOR  YES  OR  "O"  FOR  NO  :   '); 
READLN(CHSIM) ; 
WRITELN; 
IF  CHSIM  =  1  THEN 

BEGIN 

WRITELN CHOW  MANY  NUMBERS  DO  YOU  WANT  TO  USE  AS  BASIS 
FOR  THE  SIMULATION?'); 
READLN(SIMMAX) ; 

END; 
WRITELN; WRITELN; 
END; 

PROCEDURE  OUTPUT; 
BEGIN 

WRITE ('FOLLOWING  ',K:3,'  LAMBDAS  WERE  ENTERED  :'); 

FOR  I:=l  TO  K  DO 
WRITE(LCI3:5:2) ; 

WRITELN; 
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WRITE ('  #  OF  EACH  LAMBDA  :'); 

FOR  I:=l  TO  K  DO 

WRITE<NINITCI3:5)  ; 
WRITELN; WRITELN; WRITELN; 
end; 

PROCEDURE  INITIALIZE; 
BEGIN 
KK : = 1 ; 

WHILE  KK<=K  DO 
BEGIN 

ACKK, 13: =1.0; 
NICKK3:=1; 
JJ:=1; 

WHILE  JJ<=K  DO 
BEGIN 

IF  JJOKK  THEN 
BEGIN 

ACKK, 13:=ACKK, 1D*LCJJ3/ (LC J J3-LCKKD > ; 
END; 
JJ:=JJ+1; 
END; 
KK:=KK+1; 
END; 
KK:=K; 
IE:=1; 
END; 

PROCEDURE  SUMMATION; 
BEGIN 

SUM: =0.0; 
FOR  I : = 1  TO  K  DO 
BEGIN 

IF  IOIE  THEN 
BEGIN 

SUM:=SUM+ACI, 1 3*LE 1 1/ (LC I 3-LE IE3 >  ; 
IF  NICI3>=2  THEN 
BEGIN 

FACT:=1; 
NKK:=NICI3; 
FOR  II:=2  TO  NKK  DO 
BEGIN 

FACT:=FACT* (II-l) ; 
IF  LCI3>LCIE3  THEN 

SUM: =SUM+AC 1 , 1 1 D *FACT*LE IE3/EXP ( 1 1 *LN ( 
LCI3-LLIE3) ) 
ELSE 
BEGIN 

IF  (II  MOD  2)=0  THEN 

SUM: =SUM+AC I , 1 1  3 tFACTHLC IE3/EXP ( 
II*LN(ABS(LEI3-LCIE3>  > ) 
ELSE 
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SUM:=SUM+ACI, 1 1 :*FACT*LC IE3/ (-EXP( 
II*LN(ABS(LCI3-LCIE3) ) ) ) 
END; 
END; 
END; 
END; 
END; 

AC  IE, 13:=SUM+ACIE, 13; 
END; 

PROCEDURE    SUB1; 
BEGIN 
I:=l; 

WHILE    K=K    DO 
BEGIN 

IF    IOIE    THEN 
BEGIN 

ACI,NIEI33:=AEI,NICI33*LEIE3/ CLEIE3-LE 1 3 )  ; 
IF    NIEI301    THEN 
BEGIN 
J:=l; 

NNI:=NICI3-J; 
WHILE    J<=(Nini3-l)     DO 
BEGIN 

ACI,NNI3:  =  (LCIE3*ACI,NNi:-NNI*ACI,NNI-M3)  / 
(LCIED-LCI3) ; 
J:=J+1; 
NNI;=Ni:n-J; 
END; 
END; 
END; 
X:*I+1; 
END; 
END; 

PROCEDURE  ONEATATIME; 
BEGIN 

IF  NINITCIE3=NICIE3  THEN 

IE:=IE+1 
ELSE 
BEGIN 

NICIE3:=NICIE3+1; 

J:=0; 

REPEAT 

NN:=NICIED-J; 

ACIE,NN3:=LCIE3*ACIE,NN-13/ (NN-1) ; 
J:=J+1; 
UNTIL  NN=2; 
SUMMATION; 
SUB1; 
END; 
END; 
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PROCEDURE  RAND; 
VAR 

1 , 1 1 , 1 1 1 , COUNT ,  J :  I NTEGER; 
R,TSIM:REAL; 
BEGIN 

III:=0; 
COUNT: =0; 
REPEAT 

TSIM:=0.0; 
II:=1; 

FOR  J:=l  TO  K  DO 
BEGIN 

FOR  I:=l  TO  NINITCJ3  DO 
BEGIN 

R:=RANDOM/MAXINT; 
TSIM:=TSIM-M-(LN(R)  /LCII3)  )  ; 
END; 
1 1 :  =  1 1  -»- 1 ; 
END; 
IF  TSIM>TCITD  THEN 

COUNT: =C0UNT+1; 
III:=III+1; 
UNTIL  III=SIMMAX; 
RELSIMC IT3:=C0UNT/SIMMAX; 
END; 

PROCEDURE  SUB2; 
BEGIN 

RELI:=0.0; 

IF  TCITD=0  THEN 

REL I : = 1 . O 
ELSE 
BEGIN 

FOR  I:=l  TO  K  DO 
BEGIN 

SUM: =0.0; 

NKK:=NICI3; 

FOR  J:=l  TO  NKK  DO 

SUM:=SUM+ACI, JD*EXP<  ( J-l  )  *LN  (TC  ITD  )  )  ; 
PR0CID:=SUM*EXP(-LCI3*TCIT:) ; 
END; 
FOR  I:=l  TO  K  DO 

REL I : =REL I +PROL 1 1 ; 
END; 
IF  RELIM.O  THEN 

RELCITD:=1.0 
ELSE 

RELniT3:=RELI; 
IF  CHSIM=1  THEN 
RAND; 
END; 
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PROCEDURE  YAXIS; 

BEBIN 

IF  (X=0)  THEN 

(<X=0)  AND 

(X=-3>  AND 

(X=-3)  AND 

<X=-2)  AND 

(X=-l)  AND 

(X=-l)  AND 

(X=-l)  AND 

(X=-l)  AND 

(X=-l)  AND 

(X=-l)  AND 

(X=-l>  AND 

(X=-l)  AND 

(X=-l)  AND 

(X=-l)  AND 


IF 

IF 
IF 
IF 
IF 

IF 
IF 
IF 
IF 
IF 
IF 
IF 
IF 
IF 
END; 


THEN  PRINTVAR:=ASTERISK; 


PRINTVAR: ='0' ; 


PRINTVAR:=LINE; 

(Y=20))  AND  (TCI  3=0) 

<Y  MOD  2=0)  THEN 

(Y=20)  THEN  PRINTVAR:=' 1* ; 

(Y  MOD  2=0)  THEN  PRINTVAR: =  ' . ■ ; 

( (Y=20)  OR  (Y=0))  THEN  PRINTVAR: =  ' 0' ; 

(Y=18)  THEN  PRINTVAR:='9' 

<Y=16)  THEN  PRINTVAR:='8' 

(Y=14)  THEN  PRINTVAR: ='7' 

(Y=12>  THEN  PRINTVAR:='6' 

(Y=10)  THEN  PRINTVAR:='5' 

(Y=8)  THEN  PRINTVAR: =' 4' ; 


(Y=6)  THEN 
(Y=4)  THEN 
(Y=2)  THEN 


PRINTVAR 
PRINTVAR 
PRINTVAR 


=  '  1* 


PROCEDURE  GRAPHICS; 
BEGIN 

COUNT: =1; 

FOR  Y:=20  DOWNTO  1  DO 
BEGIN 

FOR  X:=-3  TO  76  DO 
BEGIN 

IF  COUNT<=TMAX  THEN 
BEGIN 

VAL : =ROUND  <  20*REL E  COUNT 1 ) ; 
IF  (VAL=Y)  AND  (X=TCC0UNT3)  THEN 
BEGIN 

PR  I NT VAR : = ASTER I SK ; 
COUNT: =C0UNT+1; 
END 
ELSE 

PRINTVAR:=BLANK; 
END 
ELSE 

PR I NTV AR : =BLANK ; 
YAXIS; 

WRITE (PRINTVAR) ; 
END; 
END; 
Y:=0; 

FOR  X:=-3  TO  76  DO 
BEGIN 

IF  (X  MOD  5  =0)  THEN  PRINTVAR: =L I NE; 
IF  COUNT<=TMAX  THEN 
BEGIN 

VAL : =ROUND ( 20*REL  L COUNT 1 )  ; 

IF  (VAL=Y)  AND  ( X=TCC0UNT3 )  THEN 
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BEGIN 

PR I NT  VAR : = ASTER I SK ; 
COUNT: =C0UNT+1; 
END 
ELSE 

IF  (X  MOD  5  <>  0)  THEN 
PRINTVAR:=DASH; 
END 
ELSE 

IF  (X  MOD  5  <>  O)  THEN 
PRINTVAR:=DASH; 
IF  (X=-3)  THEN  PRINTVAR^'O'  ; 
IF  (X=-2)  THEN  PRINTVAR:=' . ' ; 
IF  (X  =  -l)  THEN  PRINTVAR:='0'  ; 
WRITE (PR I NT VAR) ; 
END; 

O    5    10    15    20 
60'  >  ; 

TIME  — >' ) ; 


WRITE (' 

50    55 
WRITELNC 
END; 


3U 
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BEGIN 
INPUT; 
OUTPUT; 

FOR  IT:=1  TO  TMAX  DO 
BEGIN 

INITIALIZE; 
REPEAT; 

ONEATATIME; 
UNTIL  IE=K+1; 
SUB2; 
IF  CHSIM=1  THEN 

WRITELNCRELIABILITY  (  T  >'  ,  TC  111 :  3,  '  )  =  '  ,  RELC  ITD  :  S:  6, 
(SIM  =  ' ,RELSIMCIT3:5:3, ' ) * ) 
ELSE 

WRITELNCRELIABILITY  (  T  >',TEIT3:3,'  )  =  '  ,  RELE  IT3  :  9:  6)  ; 
END; 
WRITELN;WRITELN; 
ATTENTION; 

WRITELN C DO  YOU  ALSO  WANT  A  GRAPHICAL  PRESENTATION  OF 
THE  DATA?' ) ; 

WRITE ('HIT  "l"  FOR  YES  OR  "0"  FOR  NO  :   '); 
READLN(CHGR)  ; 
WRITELN; 
WRITELN; 
IF  CHGR=1  THEN 

GRAPHICS; 
DECIDE2; 
END. 
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APPENDIX  B 


PROGRAM  MIXED_CONVOLUTION; 

USES  TRANSCEND, APPLESTUFF; 
rnKiCT 

ASTERISK^  *'  ; 

BLANK    =  '  '  ; 

DASH     ='-'; 

LINE     ='  !  ' ; 

MAX      =10; 

TYPE 

ARR  =  PACKED  ARRAY C 1 .  . MAX  3 

OF  REAL; 

ARI  =  PACKED  ARRAY  LI- .MAX] 

OF  INTEGER; 

DARR  =  PACKED  ARRAY C 1 .. MAX  . 

,  1. .MAX  3  OF  REAL; 

VAR 

AR:PACKED  ARRAYCO.  . MAX, O. . MAX , O. . MAX  3  OF  INTEGER; 

A: PACKED  ARRAYC 1 . . MAX ,  1 . . MAX  3  OF  REAL; 

SUMREL: PACKED  ARRAYCO. . MAX  3  OF  REAL; 

NrPACKED  ARRAYCO. . MAX3  OF  INTEGER; 

NINIT,NI,T,ZA, INP:ARI; 

PROB , RELTOT , LAMB : DARR ; 

L,PRC,REL,P:ARR; 

I, J. K, BR, BRANCHES, COUNT, KK, IL,KKK, JJ, IE: INTEGER; 

I T . I  I , FA , NN , NKK , NN I , CHGR , TMAX , Y , X , VAL : I NTEGER ; 
SUM, RELI , LAM, DEN: REAL; 
PR I NTVAR , CH: CHAR ; 

<*  L  =  ARRAY  OF  ALL  LAMBDA  *> 

U  PROCID  =  PROBABILITY  OF  EACH  BRANCH  *> 

(*  T  =  ARRAY  OF  ALL  TIMES  AT  WHICH  THE  RELIABILITY  WILL 
BE  COMPUTED  *) 

<*  BRANCHES  =  NUMBER  OF  BRANCHES  IN  SYSTEM  *) 

PROCEDURE  INPUT; 
BEGIN 

WRITELNCHOW  MANY  DISSIMILAR  LAMBDAS  ARE  IN  THE  SYSTEM?   '); 

READLN(K) ; 

WRITELNCWRITE  DOWN  ALL  THE  LAMBDAS  WITH  AT  LEAST  ONE  BLANK 

IN  BETWEEN: ' ) ; 

WR I TELNC  (NUMBERS  <  1  NEED  A  LEADING  ZERO!)'); 

FOR  I:=l  TO  K  DO 
READ (LCI  3) ; 

WR I TELNC HOW  MANY  TIME  POINTS  DO  YOU  WANT  TO  BE  COMPUTED? 

(  <  =  20  )  :  » ) ; 

READLN(TMAX) ; 

WR I TELNC AT  WHAT  TIME  DO  YOU  WANT  THE  SURVIVAL  FUNCTION  TO 

BE  EVALUATED?' ) ; 

WRITELNC  (INTEGER  TIMES  BETWEEN  0-60  IN  THE  SAME  UNITS  AS 

THE  LAMBDAS) ' ) ; 

FOR  I:=l  TO  TMAX  DO 
READ (TCI  3) ; 
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WRITELN;WRITELN; 

END; 

PROCEDURE  ATTENTION; 

VAR 

DURATION 

, Is  INTEGER; 

PITCH: ARRAYC1. . 103 

OF 

INTEGER; 

BEGIN 

PITCHC13; 

:«  15; 

PITCHC23 

:=  17; 

PITCHC33. 

:=  19; 

PITCHC43 

:=  20; 

PITCHC53. 

DURATION' 

:  =  30; 

FOR  I:=l 

TO  5  DO 

NOTE (PITCHC I  3, DURATION)  ; 

END; 

PROCEDURE  OUTPUT 1; 
BEGIN 

WRITE ('FOLLOWING  ',K:3,'  LAMBDAS  WERE  ENTERED  :'); 

FOR  I:=l  TO  K  DO 
WRITE(LCI3:5:2) ; 

WRITELN; 
END; 

PROCEDURE  0UTPUT2; 
BEGIN 

FOR  I:=l  TO  BRANCHES  DO 
BEGIN 

FOR  J:=l  TO  MAX  DO 
BEGIN 

FOR  K:=l  TO  MAX  DO 
BEGIN 

IF  (ARC  I,  J,K3=0>  AND  (KOI)  THEN 
ARCI, J,K3:=-1; 
END; 
END; 
END; 
END; 

PROCEDURE  GETDATA; 
BEGIN 

FOR  I:=0  TO  MAX  DO 
FOR  J:=0  TO  MAX  DO 
FOR  K:=0  TO  MAX  DO 
ARCI, J,K3:=-1; 
NC03:=1; 
WRITELN; 

WRITELN  CHOW  MANY  BRANCHES?'); 
READLN (BRANCHES) ; 
I:=0; 
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REPEAT 
I:=I+1; 

WRITELN(I:3,,  •  .  BRANCH:  '  >; 
J:=0; 
REPEAT 
J:=J+1; 

WRITELN(J:3, ' .  ELEMENT:  '); 
K:=0; 
REPEAT 
K:=K+1; 
READ<2ACKD) ; 
ARCI, J,K3:=ZACK:; 
NEJ3:=K; 
UNTIL  ARCI,J,KD=0; 
UNTIL  ( (AREI, J-l,NCJ-133-0)  AND  ( ARE I , J , NE J ] 1=0) ) 
AND  (ZAE13=0>; 
UNTIL  I=BRANCHES; 
0UTPUT2; 
END; 

PROCEDURE  COMPPROB; 
BEGIN 

POR  I:=l  TO  MAX  DO 
FOR  J:=l  TO  MAX  DO 
PROBE  I, J3:=-l; 
I:=l: 
J:=2; 

WHILE  K=BRANCHES  DO 
BEGIN 
DEN:=0; 
K:=2; 
COUNT: =1; 

WHILE  AREI,J,K3>0  DO 
BEGIN 

WHILE  AREI,J,K3>0  DO 
BEGIN 

DEN: =DEN+LE ARE  I , J,K33; 
K:=K+1; 
END; 
PROBEI.COUNT3:=LEAREI, J, 13  3/DEN; 
COUNT :=COUNT+l; 
J:=J+2; 
K:=2; 
DEN:=0; 
END; 
I:=I+1; 

END; 
FOR  I:=l  TO  BRANCHES  DO 
BEGIN 

J:=l; 

PEI3:=1; 
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WHILE  probci,j:>o  do 

BEGIN 

PC  1 1 : =PC I 3XPROBC I , J 1 ; 
J:=J+1; 
END; 
END; 
WRITELN; 
END; 

PROCEDURE  COMPLAMB; 
BEGIN 

POR  I:=l  TO  MAX  DO 
FOR  J:=l  TO  MAX  DO 
LAMBCI, JJ:=-1; 
I:=l; 
J:=l; 

WHILE  K=BRANCHES  DO 
BEGIN 

LAM:=0; 
K :  =  1 ; 
COUNT: =1; 

WHILE  ARCI,J,K3>=0  DO 
BEGIN 

WHILE  AREI,J,K3>0  DO 
BEGIN 

LAM :  =LAM-t-L  t  AR  C  I ,  J  ,  K 1 1 ; 
K:=K+1; 
END; 
LAMBC I , COUNT  3 : =LAM; 
COUNT:  =CQUNT-t-l; 
J:=J+2; 
K:=l; 
LAM:=0; 
END; 
I:=I+1; 
J:=l; 
END; 
END; 

PROCEDURE  SORT; 
VAR 

TEMP: REAL; 
SORTED: BOOLEAN; 
BEGIN 
REPEAT 
I:=l; 

SORTED: =TRUE; 
WHILE  LEI+13>0  DO 
BEGIN 

IF  LEI3X.CI  +  13  THEN 
BEGIN 

TEMP:=LCI3; 
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LEI3:=LEI+13; 

LEI+13:=TEMP; 

SORTED: =FALSE; 

END; 

I:=I-H; 

END; 

UNTIL  SORTED=TRUE; 

LEI+13:=-1; 

LCI+2]:=-l; 

END; 

PROCEDURE  INITIALIZE; 

BEBIN 

KK : = 1 ; 

WHILE  KK<=K  DO 

BEGIN 

AEKK, 13:=1.0; 

NIEKK3:=1; 

JJ:=1; 

WHILE  JJ<=K  DO 

BEGIN 

IF  JJOKK  THEN 

BEGIN 

AEKK. 13:=AEKK, 13*LEJJ3/ CLE  J J 3-LEKK3 ) ; 

END; 

JJ:=JJ+1; 

END; 

KK:=KK+1; 

END; 

KK:=K; 

IE:=1; 

END; 

PROCEDURE  SUMMATION; 
BEGIN 

SUM: =0.0; 
FOR  I:=l  TO  K  DO 
BEGIN 

IF  IOIE  THEN 
BEGIN 

SUM:=SUM+AEI, 13*LEI3/(LEI3-LEIE3>; 
IF  NIEI3>=2  THEN 
BEGIN 
FA:=1; 
NKK:=NIEI3; 
FOR  I  I : =2  TO  NKK  DO 
BEGIN 

FA:=FAt (II-l) ; 

IF  LEI3>LEIE3  THEN 

SUM: =SUM+AE 1,113 <FA*LE IE 3 /EXP ( 1 1 *LN ( 
LEI3-LEIE3) > 
ELSE 
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BEGIN 

IF  (II  MOD  2)=0  THEN 

SUM:=SUM+ACI, 1 1  3 *FA*LC IE3/EXP ( 
II*LN(ABS(LCI3-LCIE3>  > ) 
ELSE 

SUM:=SUM+ACI, 1 1 3*FA»LC IE3 / <-EXP( 
II*LN(ABS(LCI3-LtIED> ) ) ) 
END; 
END; 
END; 
END; 
END; 

AC  IE, 13:=SUM+ACIE, 13; 
END; 

PROCEDURE  SUB1; 
BEGIN 
I:=l; 

WHILE  K=K  DO 
BEGIN 

IF  IOIE  THEN 
BEGIN 

ACI.NICI33:=A[:i,NICI33*L[:iE3/(LCIE3-L:i3); 
IF  NICI301  THEN 
BEGIN 
J:=l; 

NNI:=NICI3-J; 
WHILE  J<«<NICI3-1)  DO 
BEGIN 

ACI,NNI3:  =  (LCIE3*ACI,NNI3-NNI*Ai:i!,NNI  +  13)  / 
(LCIE3-LCI3) ; 
J:=J+1; 
NNI:=NICI3-J; 
END; 
END; 
END; 
I : « I + 1 ; 
END; 
END; 

PROCEDURE  ONEATATIME; 
BEGIN 

IF  NINITCIE3=NItIE3  THEN 

IE:=IE+1 
ELSE 
BEGIN 

NICIE3:=NIEIE3+1; 

J:=0; 

REPEAT 

NN:=NICIE3-J; 

ACIE,NN3:=LCIE3*ACIE,NN-13/ (NN-1) ; 
J:=J-«-l; 
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UNTIL  NN=2; 
SUMMATION; 
SUB1; 
END; 
END; 

PROCEDURE  SUB2; 
BEGIN 

RELI:=0.0; 

IF    TEIT3=0    THEN 

RELI:=1.0 
ELSE 
BEGIN 

FOR    I:=l    TO    K    DO 
BEGIN 

SUM: =0.0; 

NKK:=NICI3; 

FOR    J:=l    TO    NKK    DO 

SUM:«SUH+ACI, J3»EXP( < J-l > *LN (TE IT3 )  )  ; 
PROCI3:«SUM*EXP(-LCI3*TEITD>  ; 
END; 
FOR    I:=l    TO    K    DO 

RELI:=RELI+PR0EI3; 
END; 
IF    RELIM.O    THEN 

RELEIT3:=1.0 
ELSE 

RELEIT3:=RELI: 
END; 

PROCEDURE  PREPFORLOOP; 
BEGIN 

FOR  I:=l  TO  MAX  DO 

NINITEI3:=-1; 
KKK:=1; 
I:=l; 

IF  LE23<0  THEN 
BEGIN 

NINITE13:=1; 
K:=l; 
END; 
WHILE  LEI+13>0  DO 
BEGIN 

IF  LEI3<LEI+13  THEN 
BEGIN 

NINITEKKK3:=1; 
K:=KKK; 
KKK:=KKK+1; 
I:=I+1; 

IF  LCI+13<0  THEN 
NINITEKKK3:=1; 
K:=KKK; 


48 


IF  (X=-l)  AND  (Y=8>  THEN 

PRINTVAR:='4'  : 

IF  (X=-l>  AND  (Y=6)  THEN 

PRINTVAR: =' 3'  ; 

IF  (X=-l)  AND  (Y=4)  THEN 

PRINTVAR: =* 2' ; 

IF  (X=-l>  AND  (Y=2)  THEN 

PRINTVAR: =' 1* ; 

END; 

PROCEDURE  GRAPHICS; 

BEGIN 

COUNT: =1; 

FOR  Y:=20  DOWNTO  1  DO 

BEGIN 

FOR  X:=-3  TO  76  DO 

BEGIN 

IF  COUNT<=TMAX  THEN 

BEGIN 

VAL : =ROUND ( 20*REL  t  COUNT  3  > ; 

IF  <VAL=Y)  AND  ( X=TECOUNT: ) 

BEGIN 

PRINTVAR 

:=ASTERISK; 

COUNT :=COUNT+l; 

END 

ELSE 

PRINTVAR: =BLANK; 

END 

ELSE 

PRINTVAR: =BLANK; 

YAXIS; 

WRITE (PRINTVAR) ; 

END; 

END; 

Y:=0; 

FOR  X:=-3  TO  76  DO 

THEN 


BEGIN 

IF  (X  MOD  5  =0)  THEN  PRINTVAR: =LINE; 
IF  COUNT<=TMAX  THEN 
BEGIN 

VAL: =R0UND(20*RELC COUNT U) ; 
IF  (VAL=Y)  AND  ( X=TC COUNT 1 )     THEN 
BEGIN 

PR I NT VAR : =ASTER I SK ; 
COUNT: =C0UNT+1; 
END 
ELSE 

IF  (X  MOD  5  <>  O)  THEN 
PRINTVAR:=DASH; 
END 
ELSE 

IF  (X  MOD  5  <>  0)  THEN 
PRINTVAR: =DASH; 
IF  (X=-3>  THEN  PRINTVAR:='0' ; 
IF  (X=-2)  THEN  PRINTVAR: =' . ' ; 
IF  (X=-l)  THEN  PRINTVAR^'O'  ; 
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END 
ELSE 
BEGIN 
IIs-I; 

COUNT: =1; 

WHILE  LCIX3«LCXI+13  DO 

BEGIN 

COUNT: =COUNT+l; 

N I N I T  C  KKK 3 : =COUNT ; 

FOR  IL:=II  TO  9  DO 

LCILD:=LEIL+i:i; 
LriiAX3:=-l; 
K:=KKK; 
END; 
KKK:=KKK+1 ; 
I: =11+1; 

IF  <(<C0UNT>1)  AND  (LEI+13<0)  AND  (LEI3>0)))  THEN 
BEGIN 

NINITCKKK3:=1; 
K:=KKK; 
END; 
END; 
END; 
END; 

PROCEDURE  LOOP; 
BEGIN 

FOR  IT: =1  TO  TMAX  DC 
BEGIN 

INITIALIZE; 
REPEAT 

ONEATATIME; 
UNTIL  IE=K+1; 
SUB2; 

REL C I T D : =REL C I T 3 * P [ BR 3  ; 
RELTOTCBR. IT3 : =RELL ITD ; 
END; 
END; 

PROCEDURE  YAXIS; 
BEGIN 

IF  (X=0)  THEN  PRINTVAR:=LINE; 

IF  <<X=0)  AND  <Y=20))  AND  (TC1D=0)  THEN  PRINTVAR: =ASTER I SK; 

IF  (X=-3)  AND  (Y  MOD  2=0)  THEN  PRINTVAR: =  ' O7 ; 

IF  (X=-3)  AND  (Y=20)  THEN  PRINTVAR: =  '  1  * ; 

IF  (X=-2>  AND  (Y  MOD  2=0)  THEN  PRINTVAR: ='.' ; 

IF  (X=-l)  AND  (<Y=20)  OR  (Y=0) )  THEN  PRINTVAR: =' O' ; 

IF  (X=-l)  AND  (Y=18)  THEN  PRINTVAR: =' 9' 

IF  (X=-l)  AND  (Y=16)  THEN  PRINTVAR: =' 8' 

IF  <X=-1)  AND  (Y=14)  THEN  PRINTVAR: =' 7' 

IF  (X=-l)  AND  (Y=12)  THEN  PRINTVAR: =' 6' 

IF  (X=-l)  AND  (Y=10)  THEN  PRINTVAR: =' 5' 
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WRITE(PRINTVAR)  ; 
END; 
WRITE  ('    O    5    10    15    20    25    30    35    40    45 
50    55    60'); 
WRITELNC     TIME  — >'  >  ; 
END; 

PROCEDURE  DECIDE; 
BEGIN 

WRITELN;WRITELN; 

WRITELNCDO  YOU  WANT  TO  CONTINUE  WITH  MIXED  CONVOLUTIONS? 

TYPE  "Un  :  '  >  ; 

WRITELNCOR  DO  YOU  WANT  TO  COMPUTE  CONVOLUTIONS? 

TYPE  "XP1"  : ' ) ; 
END; 

BEGIN 
INPUT; 
OUTPUT 1; 
GETDATA; 
0UTPUT2; 
COMPPROB; 
COMPLAMB; 
FOR  I:=l  TO  MAX  DO 

LCI3:=-1; 
FOR  BR:=1  TO  BRANCHES  DO 
BEGIN 
J  :  =  1 ; 

WHILE  LAMBC BR, J  3 >=0  DO 
BEGIN 

IF  LAMB C BR, 1 3 >0  THEN 
BEGIN 

LCJJ:=LAMBIIBR,  J3; 
J:=J+1; 
END 
ELSE 
BEGIN 

LLJ3:=LAMBi:BR,  J  +  13; 
J:=J+1; 
END; 
END; 
SORT; 

PREPFORLOOP; 
LOOP; 
END; 

FOR  I:=l  TO  TMAX  DO 
BEGIN 

SUMRELCI3:=0; 

FOR  BR:=1  TO  BRANCHES  DO 

SUMRELC I  3 : =RELTOTCBR, I 3+SUMRELt I  3  ; 
WRITELN( 'RELIABILITY  (  T  >  ' ,TCI3:3, '  )  =  ', 
SUMRELEI3:8:6) ; 
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END; 

WRITELN; 
ATTENTION; 

WRITELN ('DO  YOU  ALSO  WANT  A  GRAPHICAL  PRESENTATION 
OF  THE  DATA?' ) ; 

WRITE  CHIT  "1"  FOR  YES  OR  "O"  FOR  NO  :   '); 
READLN(CHGR) ; 
WRITELN; 
WRITELN; 
IF  CHGR=1  THEN 

GRAPHICS; 
DECIDE; 
END. 
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